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Abstract 

We consider a modified form of gravity in which the action contains 
a power a of the scalar curvature. It is shown how the presence of a 
bulk viscosity in a spatially flat universe may drive the cosmic fluid 
into the phantom region (w < —1) and thus into a Big Rip singularity, 
even if it lies in the quintessence region (w > —1) in the non- viscous 
case. The condition for this to occur is that the bulk viscosity contains 
the power (2a — 1) of the scalar expansion. Two specific examples are 
discussed in detail. The present paper is a generalization of the recent 
investigation dealing with barrier crossing in Einstein's gravity: I. 
Brevik and O. Gorbunova, Gen. Relativ. Grav. 37, 2039 (2005). 



1 Introduction 

Modified versions of Einstein's gravity continue to attract interest. The mo- 
tivation for this kind of generalization has its root in the well known obser- 
vations of redshifts of type la supernovae [H El 0] , as well as the anisotropy 
of the microwave background |H El Ej- The data may be explained by dark 
energy which in turn may result from a cosmological constant, a cosmic 
fluid with a complicated equation of state, a scalar field with quintessence 
or phantom-like behaviour, or perhaps by some other kind of theory. For a 
review of the developments up to 2003, see Ref. [7j; some more recent papers 
are Refs. jHUHlEl- A forthcoming extensive review is [TTj . 

The simplest way of dealing with the expansion of the universe mathe- 
matically, is to allow for a cosmological constant. It corresponds to putting 
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the parameter w equal to -1 in the equation of state for the cosmic fluid, 



P = (j-i)p = wp. (i) 

The fluid is then the same as the extreme tensile stress vacuum "fluid" in 
the de Sitter universe. More general types of fluid can be envisaged: thus 
the region — 1 < w < — 1/3 corresponds to a quintessence fluid, whereas the 
region w < — 1 corresponds to the so-called phantom fluid. It should be noted 
that in both these cases p + 3p < 0. Thus, it follows from the Friedmann 
equation 

a 4nG , „ . A lnS 
(P + 3p) + o (2) 



a 3 3 

that the curve for the scale factor a is always concave upwards, when drawn 
as function of cosmic time t (it is here assumed that A > 0). The inequality 
p + 3p < breaks the so-called strong energy condition. 

Astrophysical data make it quite possible that the value of w is somewhat 
less than -1. It was not understood until rather recently that the universe 
will then have the strange property of passing into a Big Rip singularity 
in the future [121 HH1 El- In turn, this can lead to bizarre consequences 
such as negative entropy [15J. As shown in a recent paper of Nojiri et al. 
[To] one can now actually classify as many as four different types of the Big 
Rip phenomenon. It should also be noted that quantum effects which are 
important near Big Rip may actually act against its occurrence [17] . 

After providing these introductory remarks we have now come to the 
theme of the present work. First, we shall endow the cosmic fluid with 
a bulk viscosity (. Such a viscosity is obviously compatible with spatial 
isotropy of the universe. In a recent paper JH] we showed how the physically 
natural assumption of letting ( be proportional to the scalar expansion in a 
spatially flat FRW universe can drive the fluid into the phantom region even 
if it starts from the quintessence region in the non- viscous case. Next, we 
shall investigate if this nice property of the fluid still persists if we replace 
Einstein's gravity with a more complicated theory. Specifically, we shall 
assume the following modified gravity model: 



16ttG 



J d*x^(f R a + L m ). (3) 



Here fo and a are constants, and L m is the matter Lagrangian. This is 
the model recently studied by Abdalla et al. I24j . The dimension of fo is 



cm 2(<*-i)_ 
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Our question thus is: does there exist a crossing of the w = —1 barrier, 
from w > — 1 to w < — 1, when the field equations are derived from the 
action Q? The answer turns out to be yes, if ( satisfies the condition given 
by Eq. (jHHj) below. This is actually a quite natural generalization of the 
condition found earlier for the case of Einstein's gravity [THj . 

The present paper is a continuation of a recent study of viscous cosmology 
in modified gravity [251. Cf. also the very recent paper in on cosmological 
models with viscous media. 

We close this section by mentioning the following points. As discussed 
in [20] and [2Ij, a dark fluid with a time dependent bulk viscosity can be 
considered as a fluid having an inhomogeneous equation of state. Also, the 
possibility of crossing the w = —1 barrier was considered in these references. 
Another point is that string/M theory may predict a modified gravity action 
with a negative power of R [22] ■ A consistent modified gravity theory which 
describes not only the cosmic acceleration but may comply with the Solar 
System tests was formulated in [23] . 

2 The fundamental formalism 

We assume the spatially flat FRW metric, 

ds 2 = -dt 2 + a 2 (t) dx 2 , (4) 

and put A = 0. The four- velocity of the cosmic fluid is = (U ,!/ 1 ). In 
comoving coordinates, U° = 1, U % = 0. In terms of the projection tensor 
h^,u — 9fj,u + U fl U u we can write the energy-momentum tensor as 

T lw = pU li U v + (p-C9)h IUf , (5) 

assuming constant temperature as well as vanishing shear viscosity in the 
fluid. Here the scalar expansion is 6"^ = 3d /a = 3H, where H is the Hubble 
parameter. The effective pressure is defined as 

p = p-3H(. (6) 

From variation of the action (JHJ) we obtain the equations of motion [23J [25] 

- \fo 9^R a + af Q R^R a ~ l - af V^.R^ 1 

+af g flu V 2 R a ~ 1 = 8nGT, u , (7) 
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where T^ v corresponds to the term L m in (J3J). The notation is such that the 
values a = 1, f = 1 correspond to Einstein's gravity. 

We first consider the (OO)-component of this equation. We calculate 

Roo = -~, R = 6('- + ^]=6(H + 2H 2 ), (8) 
a \a or J 

as well as the second order derivatives of As T 00 = p we get after some 

algebra 

^f R a - 3af {H + H 2 )R a ~ l + 3a{a - l)foHR a ~ 2 R = 8nGp (9) 

(note that there is a printing error in a sign in the corresponding (00)- 
equation (16) in [25]). 

Similarly, we consider the (rr )-component of Eq. ((7j). Since in a coordi- 
nate basis 

R rr = (H + 3H 2 )a 2 , T rr = pg rr = pa 2 , (10) 

we get 

\kR a ~ af (H + 3H 2 )R a ~ l + a(a - l)f [2HR a - 2 R 

+ (a - 2)R a ~ 3 R 2 + R a ~ 2 R] = -MGp. (11) 

We consider next the local conservation equation for energy. An important 
property of Eq. (JZJ) is the following: taking the covariant divergence of the ex- 
pression on the left hand side one obtains zero. This was shown explicitly by 
Koivisto [213 • Consequently, we obtain the energy- momentum conservation 
equation in conventional form: 

VT F = 0. (12) 

Energy-momentum conservation thus occurs as a consequence of the field 
equations, just as in Einstein's theory. This property strongly supports the 
consistency of the modified gravity theory. The conservation equation for 
energy is now obtained by contracting Eq. (fT2"j) with 

p+(p + p)SH = %H 2 , (13) 

which can alternatively be written as 

p = -3-fpH + %H 2 . (14) 
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Let us next search for a differential equation for H = H(t), wherein 
{fojOijj} are given constant input parameters, and where the bulk viscosity 
C(t) is, to begin with, taken to be an arbitrary function of time. The natural 
way of accomplishing this task is to differentiate the left hand side of Eq. Q 
with respect to time, thereafter insert for p on the right hand side the ex- 
pression in Eq. ()14j). and finally insert for p again from Eq. (9). We obtain 
in this way 

\lf HR a + 3af H[2H - 3 7 (# + H 2 )^ 1 

+3a(a-l)foH[(3^-l)HR+R}R a ~ 2 +3a(a-l)(a-2)foHR 2 R a - 3 = 72nG(H 2 . 

(15) 

As R = 6(H + 2H 2 ), this can be regarded as a nonlinear differential equation 
for H(t). As the equation is complicated, its mathematical structure and 
physical meaning are best discussed in terms of examples. In the next section 
we will discuss the simplest alternative, /o = 1, a = 1. 

3 The case /o = 1, a = 1 

This case is Einstein's gravity. From Eqs. 0, (jlljl and (|15|) we obtain 
respectively 

3H 2 = 8nGp, (16) 

2H + 3H 2 = -8nGp, (17) 

2H + 3^H 2 = 24ttG(H, (18) 

in accordance with known results (cf., for instance, Ref. f2~7\). 

Of particular interest, as shown in [TH], is to take ( to be proportional to 
the scalar expansion through a constant, here called T\. Thus 

C = n e = 3 n H, (19) 

On thermodynamical grounds, ( has to be a positive quantity. For an ex- 
panding universe, therefore, T\ is also positive. If the following condition is 
satisfied, 

X = -7 + 247tGti > 0, (20) 

the equations of motion lead to the occurrence of a future Big Rip singularity 
in a finite time t. Thus even if we start from an initial situation where the fluid 
is non- viscous and being in the quintessence region (7 > 0), the imposition 
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of a sufficiently large bulk viscosity will will drive the fluid into the phantom 
region and thereafter inevitably into the Big Rip. 

On basis of the ansatz ()19|) the equations of motion are easily solvable as 
functions of t. We take the initial time as t = 0, and give a subscript zero 
to quantities referring to this instant. Defining the time-dependent quantity 
X = X(t) as 



X = 1 - xt^Gpo, (21) 
we can express the solutions as 



H=f-fp X-\ (22) 

P = PoX~ 2 , (23) 

p = p X~ 2 , with Po = wp , (24) 

R = 24nG( X + t)poX~\ (25) 
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showing explicitly how the singularities occur when X — > 0. 



4 The case a = 2 

We choose this the next step in complexity. The choice is moti- 

vated chiefly by mathematical convenience. It corresponds to a gravitational 
Lagrangian density proportional to R 2 in the action (J3J). 

Of main interest is now Eq. ()15j) , in order to make the appropriate deter- 
mination of the time dependence of the bulk viscosity. The equation becomes 

^f R 2 - 6/ [(3 7 - 2)H + 3 7 # 2 ]i? 

+ 6/ [(37 - 1)HR +R} = 72-kGQH. (26) 
In analogy with Eq. (|2*2*|) we make the following ansatz for H(t): 

where B is a nondimensional constant. In order for a Big Rip to occur, B 
has to be positive. Let us now take ((t) to be proportional to the cube of 
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the scalar expansion. Thus, denoting the proportionality constant by r 2 , we 
put 

( = T2 e 3 = 27t 2 H 3 . (28) 

The important point is that by inserting these expressions into Eq. (26) we 
find the time-dependent factors to drop out. What remains is an algebraic 
equation determining the value of B in terms of the given initial conditions: 

B 3 +[2 + -jj B 2 + - 7 B - —1 = (29) 

(note that the dimension of fo in this case is cm 2 '*" 1 ' = cm 2 ). 

The general analysis of this cubic equation leads to rather unwieldy ex- 
pressions. The structure of the equation is most conveniently discussed by 
means of examples. We define for convenience the quantity 

K = — - A (30) 
Jo 

4.1 Non- viscous fluid 

We consider this case for checking purposes. We now have r 2 = 0. From 
Eq. (|29*|l we obtain, when the trivial solution B = is disregarded, 

fi 2 +(2 + £)5 + ^7 = (31) 

for any value of the constant 7. Thus both for a "normal" fluid (7 > 1) and 
for a quintessence fluid (0 < 7 < 2/3) the two solutions of Eq. (|3*T|) are both 
negative (note that the product of the roots is equal to 37/2). There is thus 
no Big Rip solution in the non- viscous case when 7 > 0, as we would expect. 

4.2 Inclusion of viscosity, T2 > 0. Vacuum fluid 

Consider first the vacuum fluid case, 7 = 0. Equation (|29*|) yields 

B 3 + 2B 2 - K = 0. (32) 

If the left hand side of this expression is drawn as a function of B, we see 
that there is a local maximum at B — —4/3 and a local negative minimum 
at B — 0, irrespective of the magnitude of K. There exists thus one single 
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positive root of the equation, for all positive K. This root is caused by the 
viscosity, and leads to the future singularity. Let us assume that K (or r 2 ) 
increases from K — upwards; then we first encounter a parameter region 
in which there exist three real roots [28J. Let us assume this region in the 
following, and introduce an angle <fi in the interval < <fi < 180° such that 

/ 27 \ 

cos0 = -(l K) . (33) 



16 

The solutions can then be written in the form 

" g ~~ \~ g COS 77 5 

2 _i 4 „„„ ?<* , ir, n 



fi = <| -| + | cos § + 120" , (34) 
-! + f cos (1 + 240°). 

As an example, we choose the value K = 8/27, corresponding to <p — 120°. 
Then, B = { — 1.9196,-0.4351,0.3547}, where the last positive solution de- 
scribes the viscosity-generated Big Rip phenomenon. 

4.3 Inclusion of viscosity, t<i > 0. The fluid almost a 
vacuum fluid 

As measurements show that the parameter 7 is close to unity, it is of physical 
interest to make a perturbative expansion around 7 = 0. Let us consider the 
formalism to the first order in 7, assuming 

|7|<1- (35) 

The basic equation is Eq. (j2T?j) . as before. Drawing the left hand side of 
the equation versus B, we see that the rightmost extremal point is a local 
minimum, with coordinates (— 17, — K). Therefore, the equation must have 
one positive root. Choosing again the value K = 8/27 as an example, we can 
as before define an angle in the region < <p < 180°; it is now determined 
by 

cos0 = — - ( 1 — —7 ) . (36) 



2 V 16 
The positive root is then found to be 

leading to the viscosity-generated Big Rip in this particular case. If 7 = 0, 
the first member of Eq. (J34)) is recovered. 
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5 Remarks on the general case 



For arbitrary values of a we can no longer give complete solutions, but it 
turns out that the most important part of the problem, namely to determine 
which form of viscosity leads to a Big Rip, can easily be dealt with. First, we 
assume the same form of H = H(t) as before; cf. Eq. (|27|). where the value 
of the parameter B depends on a. This means that R = 6(H + 2H 2 ) has 
the same form as before. In analogy with Eqs. (j!9|) and (|2*Hj) we assume next 
that ( is proportional to the scalar expansion raised to the power (2a — 1). 
Denoting the proportionality constant by r a , our basic ansatz thus reads 

C = Ta 02^1 =Ta(3H) 2 a ~l_ (3g) 

Upon insertion of the expressions for H, R and £ into Eq. (j!5|) we see that 
the time-dependent factors again drop out, and we remain with the following 
equation determining the value of the constant B: 

(B + 2) a ~ 1 {9(2 - ot)7 + 3 [a + 37 + a(2a - 3)(3j - l))B 

+ 6a(a - l)(2a - 1)B 2 } - ^ (-Y nGr a = 0. (39) 

J fa V 2 J 

Because of its complexity this equation has to be analyzed in each specific 
case. Any positive root for B leads to a Big Rip singularity, as follows from 
the form fl27|) for H(t). 

If a = 2, Eq. (J39"J) reduces to our previous Eq. (J23J) . It is of interest to 
reconsider in more detail the case a — 1: then Eq. (|39j) yields the solution 

£ = --7 + 367^1. (40) 

Since in this case 

I8ttG 

#o = y— Po (41) 
according to Eq. (J22|) . it follows upon insertion that 

X = 1 - BH t = 1 - xt\/6nGp , (42) 

where x is defined by Eq. (1201) . This is agreement with our previous expres- 
sion (|21|) for X when a = 1. 
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6 Summary 



Our treatment was based on the form © for the action. Then assuming the 
cosmic fluid to possess a bulk viscosity ( varying with the scalar expansion as 
in Eq. (jSHJ), we found how the fluid could in principle be driven into the Big 
Rip singularity, even if it stayed in the quintessence region in the non- viscous 
case. This is a property following directly from Eq. (}2"Tj) for H, whenever the 
constant B is positive. In general, B is determined by Eq. We worked 

out the solutions in reasonable detail when a = 1 (Einstein's gravity), and 
also when a = 2, assuming a vacuum fluid (7 = 1 + w = 0) as well as almost 
a vacuum fluid (|7| <C 1). 

Finally, we mention that it would be of interest to incorporate the above 
formalism in the modified Gauss-Bonnet gravity interpreted as dark energy 

mum- 
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